Light hadron spectrum with Kogut-Susskind quarks* 
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We made an extensive study of the light hadron spectrum using the Wilson gauge action and Kogut-Susskind 
quarks. Using both dynamical quarks and the quenched approximation, we determine hadron masses for five and 
four gauge couplings respectively, with at least five quark masses at each coupling. In the continuum limit, we 
find a significant difference between two flavors of dynamical quarks and the quenched approximation for a range 
of quark masses. 



Calculating the light quark spectrum is a long 
term goal of QCD. Because of limitations of the 
numerical approach, one must work with quark 
masses heavier than in Nature and with a non- 
zero lattice spacing. Good control over the ex- 
trapolations in quark mass and lattice spacing re- 
quires that simulations be carried out over a wide 
range of quark mass and lattice spacing. We have 
been studying the spectrum with Kogut-Susskind 
or staggered quarks for quite some time Q . This 
year, we extended our dynamical data set and 
developed a new extrapolation method for the 
quenched case. 

In the quenched approximation, in the contin- 
uum limit at the physical value of m w /m p , we 
find m N /m p = 1.254 ± 0.018 ± 0.028 @, where 
the first error is statistical and the second is sys- 
tematic. With dynamical quarks, the double ex- 
trapolation in quark mass and lattice spacing is 
difficult. However, we can confidently extrapolate 
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to the continuum limit at intermediate values of 
the quark mass where we only need to interpo- 
late in quark mass or extrapolate slightly from 
our lightest quark mass. For dynamical quarks 
at such intermediate quark masses, mj\r/m p is 
significantly larger than it is in the quenched ap- 
proximation. We believe this is the first convinc- 
ing evidence of the effects of dynamical quarks 
on the light quark spectrum. Before having con- 
fidence in our extrapolation to the chiral limit, 
we would need to have results for lighter quarks 
at our two weakest couplings. At Edinburgh j|, 
we presented our continuum extrapolation for the 
physical value of m^/m p and found it consistent 
with the experimental value. However, at that 
time we only had four dynamical quark masses at 
the two weakest couplings. Having added a fifth 
mass in each case and choosing some different fits 
of hadron propagators, some of the extrapolated 
masses have changed by two standard deviations. 
Thus, we urge some caution not to overreact to 
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our dynamical results in the chiral limit. 

Some details of the configuration generation 
and propagator calculation may be found in 
Ref. jp. For the dynamical runs, the minimum 
values of m-n/rrip are 0.48 and 0.53 for 6/g 2 — 5.5 
and 5.6, respectively. For the stronger couplings, 
we have gone below 0.4. 

We have found the chiral extrapolation for 
the quenched simulations to be a most vexing 
problem. Quenched chiral perturbation theory 
(QxPT) indicates that there are corrections to 
the nucleon and p masses that do not occur with 
dynamical quarks. The chiral expansion is com- 
monly expressed in terms of the quark mass. For 
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the quenched case, the terms m q and m q log m q 
are the new terms. (They are due to rj' loops.) 
However, we note that due to flavor symmetry 
breaking for staggered quarks, the flavor singlet 
pion that appears in QxPT does not actually have 
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a mass proportional to m q . Thus, it is actu- 
ally more appropriate to express the chiral expan- 
sion using a term proportional to a non-Goldstone 
pion mass. We define, for fixed Ajv, 



phys 

m N = (m N + X N m n2 ) ph ^ 

mr^ + ANrnZ 



(1) 



where phys stands for the physical values and 
the other quantities are values computed at a 
given quark mass and lattice spacing. A simi- 
lar equation applies to the p. We then fit m' N 
and m' p using only terms that appear in ordi- 
nary xPT (M + am q + bm 2 plus either crr? q 2 
or cTO^logmg) for various values of Ajv and A p 
obeying < \n < A p < 0.4, which is the range 
expected from QxPT. Further details of the chi- 
ral and continuum extrapolations may be found 
in Ref. [|. 

For the dynamical quark calculations, we tried 
six fits that only contain terms appropriate to or- 
dinary xPT. The same four parameter forms we 
used above give good combined confidence level 
for the nucleon of 0.55 and 0.41. However, for the 
p the combined CL is quite poor. The poor fits 
occur for the strongest couplings. At 6/g 2 = 5.3, 
in particular, there are two points that lie off a 
smooth curve. So far, we have picked fits inde- 
pendently at each run; however, we expect that 
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Figure 1. Continuum extrapolation for m^/mp 
with the physical quark mass. 

the appropriate fitting range should be a smooth 
function of the quark mass. We may adjust our 
fitting ranges with this in mind. 

The final extrapolation is in the lattice spac- 
ing. The leading error for Kogut-Suskind quarks 
is 0(a 2 ). For the quenched case, we must consider 
different values of A^r and lambda p , and the sys- 
tematic error comes from this variation, as well as 
considering a possible quartic lattice spacing er- 
ror, or a contamination from the iti that results 
in a linear error. When we combine all these sys- 
tematic errors in quadrature, we find our final 
result m N /m p = 1.254 ± 0.018 ± 0.028. 

For the dynamical quark results, Xn = A p = 0, 
so the analysis is more straightforward. Note, 
however, that the lattice spacing changes with 
quark mass. It is therefore helpful to fix m^/m p 
(and hence the quark mass) to various values 
and then to do the continuum extrapolation sep- 
arately for each choice of m n /m p . We compare 
the dynamical case with the quenched by letting 
A at = A p — and using the same fitting function 
for N and p. (The terms specific to Qx?T are 
only large for very small quark mass.) In Figs. 1 
and 2. we compare quenched and dynamical re- 
sults at m^lrrip = 0.1753 and 0.5 with continuum 
fits that include an error quadratic in the lattice 
spacing. 
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Figure 2. Continuum extrapolation for m^/mp 
with m^/nip = 0.5. 
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Figure 3. Edinburgh plot in the continuum limit 
for N f = and 2. 



In Fig. 3, we compare the quenched and dy- 
namical results in the continuum limit. (The 
burst, slightly displaced to the left, includes both 
systematic and statistical errors 0.) The dy- 
namical results are above the quenched results 
for the entire range of m w /m p we consider. For 
■m^/nip = 0.55 there is clearly a significant dif- 
ference between the two (0.041 ± 0.007). This 
is a value where we don't have to extrapolate in 
quark mass. Also, if the quark mass were really 
this heavy, the p would not decay. We are seeing 
a real difference between the two flavor dynam- 
ical spectrum and the quenched approximation. 
Of course, we are especially interested in the chi- 
ral limit, because there we can compare with the 
real world. Recognizing that the quenched result 
for the physical value of m^/m p is closer to the 
experimental result, we are left to speculate why 
this is the case. First, we wish to again mention 
that for the dynamical case we have less confi- 
dence in the chiral limit because of the long ex- 
trapolation for the two weaker couplings. Sec- 
ond, as has always been the case with dynamical 
simulations, we take no special measures to deal 
with the potential for p decay. In lattice simula- 
tions, the threshold for decay is changed by finite 
volume effects. Perhaps the lattice suppression of 



the decay changes the p mass, or a complete anal- 
ysis of the p and two tt channels is required to cor- 
rectly determine the p mass. Third, we only have 
two dynamical quarks in our calculation. Per- 
haps the strange quark plays a larger role than 
we might have expected. 

Although much additional work remains, par- 
ticularly regarding the chiral extrapolation, we 
are starting to obtain reliable results for lattice 
QCD, even with dynamical fermions, in the con- 
tinuum limit. We see a very interesting difference 
between the Nf = 2 theory and the quenched ap- 
proximation even in a region where only chiral 
interpolation is necessary. 
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